arXiv: 1509.05182v2 [math.AP] 30 Sep 2015 


Ground State Patterns and Phase Transitions of Spin-1 
Bose-Einstein Condensates via F-Convergence Theory 

I-Liang Chern*’^’^, Chiu-Fen Chou^’^, and Tien-Tsan Shieh^’^ 

^National Center for Theoretical Science, National Taiwan University, Taipei, 106, Taiwan 
^Department of Applied Mathematics, National Chiao Tung University, Hsinchu 300, Taiwan 
^Institute of Applied Mathematical Science, National Taiwan University, Taipei, 106, Taiwan 


Abstract 

We develop an analytic theory for the ground state patterns and their phase transitions 
for spin-1 Bose-Einstein condensates on a bounded domain in the presence of a uniform 
magnetic field. Within the Thomas-Fermi approximation, these ground state patterns are 
composed of four basic states: magnetic state, nematic state, two-component state and 
three-component state, separated by interfaces. A complete phase diagram of the ground 
state patterns are found analytically with different quadratic Zeeman energy q and total 
magnetization M for both ferromagnetic and antiferromagnetic systems. Using the E- 
convergence technique, it is found that the semi-classical limits of these ground states 
minimize an energy functional which consists of interior interface energy plus a bound¬ 
ary contact energy. As a consequence, the interface between two different basic states has 
constant mean curvature, and the contact angle between the interface and the boundary 
obeys Young’s relation. 


1 Introduction 

1.1 The Gross-Pitaevskii equation for spin-1 BECs 

In the 1920s, Bose and Einstein [|9ll21 predieted a new state of matter. At very low temperature, 
very dilute Boson gases such as alkai gases tend to occupy their lowest quantum state. The de 
Broglie wave of the bosons is coherent and their lowest quantum state becomes apparent, called 
Bose-Einstein condensation (BEC). It was only recently, in 1995, BEC was first realized in lab¬ 
oratories by two groups independently, Cornell and Wieman as well as Ketterle [l2l|T8l. Through 
the mean field approximation and Hartree’s ansatz, the mean-field state of an A^-particle boson 
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gases can be described by a complex order parameter t). Its dynamics is modelled by the 
Gross-Pitaevskii equation [fT7ll25l[5T]l : 

ihOti) = “^1^) where = J H{'ip)dx, 

and the Hamiltonian of the system is written as 

H{i^) := 

Here, denotes the complex conjugate of 'i/j, and 58/5ip* is the variation of 8 with respect to 
Ip*. The function V(x) is the trap potential satisfying l/(a;) —)■ cxd as |x| —)■ cxd. The parameter m 
is the mass of the boson, and /? is the product of the particle number N and interaction strength. 
This interaction is attractive when /3 < 0 and repulsive when pJ > 0. This mean field model was 
rigorously justified as a limit of the -particle quantum system by Lieb et al. [l35l|36l for the 
ground state cases and by Erdos et al. for the dynamic cases [l22ll . 

When an optical trap is applied to confine BECs, all their hyperfine spin states are active. 
Such a BEC with an internal spin freedom is called a spinor BEC Il28ll47l[34l and was realized 
in laboratories with ^^Na atoms in 1998 ll54l l53l ISTl [8l. In the mean field theory, a physi¬ 
cal state of a spin-/ BEC is described by (2/ -f l)-components of complex order parameters 
= {ipf, i’f-i, • • •, '/o • • ■, i’-fY its dynamics is governed by a multi-component Gross- 
Pitaeviskii equation [[28l 071 . In this paper, we concern ourself with a quantum system of the 
spin-1 BEC, whose dynamics are described by a generalized Gross-Pitaevskii equation 

ihdt^ = ®, £:[vi/] = j H{^) dx, (1.1) 

where = (/>i, ipQ,ip_iY and the Hamiltonian is given in the form 


H{m) ;= Hkinp^) + + Hs{^) 

2m' ' ^ ' 2 ' ' 2 ' ' 

Here, := {'ipl,'ipQ,'ip*_i) and F = Fy, F^) is the spin-1 Pauli operator: 




/ 0 1 0 \ 
1 0 1 
Vo 1 0/ 


F = _ 


/ 0 -1 0 \ 
1 0 -1 
lo 1 0 / 


,F.= 


( 1 0 
0 0 
V 0 0 


0 \ 
0 

-1 / 


The term ^ | | ^ represents the spin-independent interaction between bosons, whereas y | fh^Ffh p 
stands for the spin-exchange interaction between bosons. The spin-independent interaction is 
attractive if c„ < 0 and repulsive if > 0. This BEC system is called ferromagnetic if c^, < 0 


2 












and antiferromagnetic if > 0. In experiments, an example of a ferromagnetic system is 
the alkali atom with c„ ~ 7.793, Cg ~ —0.0361 []3, whereas for the alkali atom ^^Na, 

Cn ~ 15.587 and Cg ~ 0.4871 [fl^ . and it is an antiferromagnetic system. 

We are interested in the ground state patterns of spin-1 BECs in the presence of an external 
uniform magnetic field. The interaction of atoms with the applied magnetic field, say Bz, 
introduces an additional energy, called the Zeeman energy: 

EzeeM = [ HzeeWdx, Hzee{^) := 

i=-i 

where Ej (B) is the Zeeman energy shift for each component under the magnetic field Bz. It is 
convenient for later discussion to write the Zeeman energy Hzee as 
1 

Hzee = Ej{B)\ij_i\^ = g(|^/>i|2 + |^/>_iH+p(|V'i|2-|'0_ip)+Eo(|V’iP + |'0on+V’-in, 
i=-i 

where 

P = \iE_i - El), q = ^{E-i + Ei- 2Eo). 

The parameters p and q are called linear and quadratic Zeeman energy, respectively, in the 
physics literature. 

Notice that the Gross-Pitaevskii equation for spin-1 BECs possess two invariants: the total 
mass 

A/'[\I>] := J\'^\‘^dx = J \'ijji\‘^ + + \^jJ-lf dx, (1.2) 

and the total magnetization 

:= J li^il^-li^.il'^dx, (1.3) 

which can be derived by direct calculation. The ground states of the system are those physical 
states of lowest energy, given fixed total mass and fixed total magnetization. This defines a 
variational problem for ground states in the presence of a uniform magnetic field: 

mm{^[^] : Af[^] = iV,7W[^] = M} . (1.4) 

where the energy takes the form 

£^[4/] = J (^Hkin + Hpot + Hn + Hg + g(|'0ip -I- |'?/’_ip)) dx + EqN + pM 

= J |^|V^|2 + V^(a;)|^|2 + ^|^|4 + ||^tF^|2 + g(|V^,|2 + |^_,|2)rfa; + ^QiV + pM. 

It is observed that the parameter p plays no role in minimization for fixed total magnetization M. 
Thus, our goal is to study the ground state patterns and their phase transitions in the parameter 
plane (g, M). 
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1.2 A brief survey of ground state problems 


There are several studies eoneemed with ground states of spin-1 BEC systems. The existenee of 
ground states for spin-1 BECs in three dimensions with |cs| < Cn was given in [|^ . The non- 
existenee result in three dimensions for Cn < 0 was given in a. Other existenee/non-existenee 
results in one dimension are given in ifTTll . 

In the ease of no applied magnetie field, it is proven that the ground state is the so-ealled 
single mode approximation (SMA) for a ferromagnetie system. That is, the ground state has the 
form ( 71 , 7 o, 7-i)'0, where 7 ^ G M"*", and -0 is a sealar field. On the other hand, the ground state 
is a so-ealled two-eomponent state (2C) for antiferromagnetie systems. The above results are 
given numerieally in [1T71[T31I . and analvtieallv in lf3^ . 

When there is an applied magnetie field, it is found that there is a phase transition from 2C 
to 3C, i.e. all three eomponents are not zeros , as q > qcr > 0 for antiferromagnetie systems. 
This phase transition phenomenon was also observed in the laboratory [l54l [33 . It was found 
numerieally in [[5911371 ITTII and proven analy tie ally in [[39ll . 

Due to the faet that the eoeffieient ^ 1, people have paid attention to the Thomas-Eermi 

approximation of ( |1.4| ) whieh simply ignores its kinetie energy. Under sueh an approximation, 
existenee results for ground states and phase transitions with |C,| « Cn have also been studied 
by many physieists, see [l59ll42ll40llTO . 

Without ignoring the kinetie energy, the problem ( |1.4| ) is a variational problem of singular 
perturbation type whieh is similar to the Van der Waals-Cahn-Hilliard gradient theory of phase 
transition for a fluid system of binary phases eonfined to a bounded domain under isothermal 
eonditions. Gurtin [l27l had eonjeetured about the asymptotie behavior of the variational model 


inf Iy u\'^ +-W (u)^ : u G J u{x) dx = m 

where IE : —)■ M+ is a double-well potential and m is the total mass of the fluid. The 

problem have been studied extensively by many authors through de Giorgi’s T-eonvergenee 
theory. The sealar ease (n = 1) is studied by [I431133ISSll . In [|24l|55l[56||, the veetorial ease 
(n > 1) is studied. With a given Diriehlet boundary eondition, a sharp interfaee limit of the 
energy funetional is eonsidered in [j^ |49l 153 . A minimal interfaee problem arising from a 
two-eomponent BEC in the regime of strong eoupling and strong segregation was studied in 
[[T]| via E-eonvergenee. They have formulated the problem in term of total density and spin 
funetions, whieh eonvert the energy into a sum of two weighted Cahn-Hilliard energy. 



1.3 Contribution of this paper 

This paper eonsiders ground state patterns and their phase transition diagram on the q-M plane 
for the ease [c^l < c„. The eontribution of this paper ineludes 
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• we find all possible ground states eonfigurations for the spin-1 BEC system in its Thomas- 
Fermi approximation and give a eomplete phase diagram on the parameter spaee (g, M); 

• a sharp interfaee limit of the BEC system is derived through de Gorgi’s E-eonvergenee. 

This paper is organized as the follows. Seetion 2 is a reformulation and simplifieation of 
the problem. Seetion 3 and 4 are devoted to the Thomas-Fermi approximations and the E- 
eonvergenee results. 


2 Formulation of the problem 


In this Seetion, the variational problem ( |1.4| ) for ground states of the spin-1 BEC system is 
formulated into a real-valued variational problem with dimensionless eoeffieients. 


Dimensional Analysis We minimize the energy funetion E ['E] under the two eonstraints (1.2), 


( |1.3| ). The eorresponding Euler-Eagrange equations read 


(/i -t- \)ipi = 

i^i>o = 

(/i - = 


■+ V {x) + q + Con 
2m 


- — -f C(x) -f Con 

2m 


■—'V‘^ + V{x)+q + Con 
2m 


+ C2{ni +no- -f C2i^*_i'ipo^ 

ijjo + 2c2(ni - n_i)V^o + C2'0i'0-iV’o 

+ C2(n_i + no- ni)'ijj_i + C2ipl'ipl. 


Here, /i and A are the two Eagrange multipliers eorresponding to the two eonstraints. We 
perform resealing: x = Lx', ip = L~^/‘^ip', then eompare the dimensions of the first equation: 


—L-^^‘^L-^V^iP' + VL-^/^iP' + qL-^l‘^i)' 4 - CnL-^n'ip'^ + CsL 
2m 


-3 


From this, we define new parameters: 


^2 _ L_j^-7/2 

2m 

V' = 

q' = qL~'^l‘^ 

~ ^nL 


-3 


c; = c,L-L 


Dropping the primes, we get the Hamiltonian defined in (2.2) below. 
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Bounded domain problem with zero potential The trap potential V in the laboratory satis¬ 
fies V(x) —)■ oo as |a;| —)■ oo, whieh leads to the exponential deeay of the ground states. This 
fact can be derived from standard elliptic PDE theory [[35l [39]l . In particular, the quadratic po¬ 
tential (which is commonly used in laboratory) is close to zero potential near the trapped center. 
Thus, it is reasonable to consider the following potential with infinite well 


V(x) 


0 if a; G 
oo if X ^ Q, 


where is a smooth bounded domain in This corresponds to the constrained variational 
problem in a bounded domain with zero Dirichlet boundary condition. 


Reduction to a real-valued problem To study the ground states, we also notice that we 
can limit ourselves to those order parameters ^jJj with constant phases. In fact, if we express 
'ipj = Uj > 0, j = 1,0, —1, then the kinetic energy is \'V'ijjj\‘^ = iVwjp -f 

which has minimal energy when VOj = 0. In this situation, the only term in the Hamiltonian 
H which involves phases is 

= y ((“i “ “-i)^ + 2'Uo(“i + “-1 + cos(A6')) 

where A9 = 9i + 9_i — 29q. The Hamiltonian Hg has a minimal value when 

A »=/0 

( TT if C5 > 0. 

The resulting Hg becomes 

Hg{u) = j ((uf - + 2uo(ui - sgn(cg)u-iy) . ( 2 . 1 ) 

Summary of the Problem To summarize the above simplifications, we shall consider the 
following constrained variational problem: 

f inf {£glu] := H,(u) dx | u e H L^(Q,Rl)} 

(Pg) < subject to the constraints 

[ A/’[u] := f^(uj + Mq + "^-i) Al[u] := f^(uj — m^j^) dx = M. 

Here, the Hamiltonian Hg(u) is 

f/,(u) := ^ \Vuj\‘^ + ^\u\'^ + ^[2ul{ui-sgn{cg)u-i)‘^ + {ul-u\)‘^] +q{ul + u\). 

i=-i 

( 2 . 2 ) 
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3 The ground states in Thomas-Fermi approximation 


The semi-classical regime is the case where e is small (for instance, choosing large L). The 
Thomas-Fermi regime is the case where e = 0. Thus, we split the Hamiltonian into 


H,{vl) = |VM,f + (3.1) 

i=-i 


A Thomas-Fermi solution is a measurable function u on H which solves the constrained varia¬ 
tional problem: 


inf {£^'ri7’[u] := f^HTF(u)dx \ u G (L^(f2, 
subject to the constraints A/'[u] = N, Ad[u] = M. 


It is expected that the Thomas-Fermi solutions are piecewise constant solutions consisting 
of one or two pure states in the following forms: 


• Nematic State (NS), if u = (0, Uq, 0) 


• Magnetic State (MS), if u = {ui, 0, 0) or (0, 0, m_i) 


• Two-component State (2C), if u = (mi, 0, m_i) 


• Three-component State (3C), if u = {ui, uq, M-i). 

Here, m* above denotes a nonzero value. We shall give a complete phase diagram of the ground 
states in this section and describe the ground state patterns in the next section. For a given 
total magnetization M, there exist two critical numbers gi, q 2 such that we have the following 
description of the phase diagram: 

• For Cs > 0,q2 < q, the Thomas-Fermi solution is a NS + MS mixed state. 


• For > 0, gi < g < g 2 , the Thomas-Fermi solution is a NS + 2C mixed state. 


• For > 0, g < gi, the Thomas-Fermi solution is a 2C pure state. 

• For < 0, g < 0, the Thomas-Fermi solution is a MS + MS mixed state. 


• For < 0, 0 < g, the Thomas-Fermi solution is a 3C pure state. 

Precisely, the notation NS + MS means that there is a measurable set (7 C H such that 

u(a;) = a-Xu{x) -f bxo\c/(a;) 

where xu and Xn\u are characteristic functions and the vectors a = (0, uq, 0) and b = (wi, 0, 0) 
or (0,0, M_i) are two constant states. 
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Normalization and Notations We may divide Htf by c„, and set a := Cg/cn and rename 
q/cn still by q. We eonsider to be a bounded set with smooth boundary. With this normaliza¬ 
tion, the Thomas-Fermi Hamiltonian beeomes 

:= + ul + u\Y + ^[2ul{ui-sgn{a) u_iY + {ul-u\Y]+q{ul+u\). (3.2) 

We denote the ratio |(7|/|f2| by r, the mass per unit volume iV/|f2| by n, and magnetization per 
unit volume M/|f2| by m, respeetively. 


It is observed that the role exehanging between ui and m_i does not ehange the form of ( |3.2| ) 
and the eonstraint of the total mass ( |1.2| ). However, it ehanges the sign of the total magnetization 
( |1.3| ). Beeause of this symmetrie property, we only need to eonsider the ease m > 0 without 
loss of generality . 

3.1 Antiferromagnetic BEC (ct > 0): ^2 < ^ implies NS + MS state 

Theorem 3.1. Suppose 0 < a < 1, m > 0. Let 

q2= (l- (n+((a + l)'/2-l)m). 

Then for q > q 2 , the global minimizer of the constrained variational problem (Pq) in the finite 
domain takes the form 

u = axc/ + bxo\c/ 

where U <Z fl is a measurable set of size 

{a + 


\U\ = 


n 


+ ((« + 1)1/2 -1) 


m 




(3.3) 


A 


a = 


(a+ 1)1/2 


,0,0), b = (0,vC4,0) 


and 


A = (n + ((a + 1)^/^ — 1 ) m) . 

Proof 1. First, we rewrite as the sum of several perfeet squares: 

2Htf = {u\ + Ug + + a{ul - + 2aul{ui - U-iY + 2q{ul + u\) 

= (1 + q)u\ + Mq + (1 + q.)N_^ + 2(1 — q)u\vS^ 

+ 2uq + u_i) + cr('Ui — uS ] P 2q{u^ + u_f) 


= ( (1 + aY^'^ul + Mq + 


a 


tU 


2Un 


(l + a)l/2--l 

(1 + a - (1 + a)^/^) ul+ (1 + 


+ uti + 2,q{u\ + u\) 


1 + 0 
1 — 0 


o 


(l + o)l/2 


u 


-1 


2auiU-i 


8 


















In the last term, the quadratic form is non-negative because 


(l -|- a — (1 -f u\ + il + a 


1 — a 

(l + a)V2 


u\ — 2auiu_i 


((a -f 1) - (a -f 1)^'^^) Ml 


a 


((« + l)-(« + l)V2)V2 


M-1 


2a + 4-2(a + l)i/2 


u\ 


> 0 , 


(a -f-1) — (a -I- 1)^2 

Here, we have used a -f-1 — (a -I- 1)^/^ > 0 and 2q; -t- 4 — 2(a -h 1)^/^ > 0 for 0 < a < 1. 


2. Because of the constraints of total mass and total magnetization, we can convert the varia¬ 
tional problem to an equivalent one by adding the terms involving j3i dx, 

/^ 2 /n {u\ — v?_i) dx and some constant dx. Namely, we can replace Htf with an¬ 
other energy density defined by 

H^(u) := HtfH + y (mi + Mq + «-i) + y K - «-i) + Y’ 

so that the minimization problem (Pq) is equivalent to 

inf / iy(u) dx subject to (1.2), (1.3). 

Jn _ _ 

We shall choose the parameters A , [32 and A so that 


(1.2),(1.3 


W{ai) = W (b) = 0, br(u) > 0 otherwise. 


(3.5) 


We organize W as 


2W= (^(a + + 


1 — a 


(l+a)i/2 

-f 2 mo j ((a -f 1) - (a -f 


M?.i - A 


a 


((«-t-l)-(« + l)V2) 


1/2 


M-1 


2a + 4-2(a + l)^/2 


2 2 I 4 


^ ^ I 

U—iUq + 


0( -f- 1 


(a -|- 1 ) — (a -f 1)^P 
+ (/?! -)- /)2 + 2q' -f 2A(o! -f 1)^^^) -f (/?! -f 274 )mq -f ^/3i — /32 -h 2g A 2A 


1 — a 

(a+ 1)1/2 


3. Now, we choose /3i and [32 to satisfy 


[3l A [32 3- 2q -\- 2y4(Q! + l)l/2 
/?! + 2^4 


= 0 
= 0 
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and the coefficient of ul 1 becomes 


/?! - /?2 + 2g + = 4A 


We define g 2 by 


Now, 


(a+ 1 ) 1/2 

g 2 := ^ - 


(a+ 1 ) 1/2 


-l] +4q. 


{a + 1 ) 1/2 


21+ = ^(a + 1)^/^M^ + "*^0 + 


1 — a 


(l + a)i/2“-i 
+ 2ul ( ((a + !)-(« + 


A 


a 


((a+ !)-(«+ 1 ) 1 / 2 )'/' 


u-i 


2a + 4 — 2(a + 1)1/^ 2 2 , ^a 4 ,^2 

Mo“-i + —^ “-1 + 4 (g - g2) 


(a + 1) - (a + 1)1/2 “0--1 ■ ^ + 1 -1 — “-!• (3-6) 

When g > g 2 , the coefficient of U-i is positive. Every term in the function W is non 
negative and the only zero of the function W{ui,uo,U-i) in is 


A 


a = 


(a+ 1 ) 1/2 


, 0,0 


and 


b = (0,\/)4,0). 


4. If u(x) = a or b for x a.e. in f2, then IV (u(a;)) dx = 0. On the other hand, since IV > 
0 , we get that if W (u(a;)) dx = 0 , then W (u(a;)) = 0 for x a.e. in O and thus u(x) = a 
or b a.e. in O. Therefore, a minimizer of the variational problem inf W (u) dx can be 
expressed as 

u = a xr/+ b (3.7) 

for some measurable set U c O. 


5. Finally, we determine the constant A and the measure |(7| by plugging (3.7) into the two 

'1 |f/| +/1(|SJ| - |(7|) = iV 


constraints (| 1 . 2 |), (11.31): 


a 


1 ) 1/2 


A 


(a+ 1 ) 1/2 


\U\ = M. 


These lead to 


A 


IC'I 


+ (/V + (a + - M) 

(n + ((a + 1 )'/' — l)m) 

(a + l)i/ 2 m 
n + ((a + 1 ) 1/2 _ 


□ 
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3.2 Antiferromagnetic BEC {a > 0): q < qi implies 2C state 


Theorem 3.2. Suppose 0 < a < 1 and m > 0. Let 


Qi = 



+ s/'o? -\- amP^ . 


(3.8) 


Then for q < qi the global minimizer of the constrained variational problem (Pq) is the constant 
state 


u = 


n + m In — m 


(3.9) 


Proof 1. Under the eonstraints of the total mass and total magnetization, the variational 
problem inf dx is equivalent to the variational problem inf IV(u(x)) dx 

under the same eonstraints, where 

W(ui,Uo, U-i) := JLtf(ui, Uo, U-i) + ^ + 

- (A + q)(ui + Mq + U-i) - aB{u\ - u\). (3.10) 


The goal is to show that lU(u) > 0 for 0 < g < gi, and the only zero of W satisfies (3.91. 
The eonstants A and B will be determined by the two eonstraints. 


2. The strategy is to introduee two parameters ki and ^2 to make the eoeffieient of Uq to be 
zero, to maximize the eoeffieient of Uq, and to make the rest to be non-negative: 

2W = («! + Mq + U-i — A)^ + a(u^ — u^_i — B^ + 2aul{ui — m_i)^ — 2gMQ 

= {ul + kiul + — AY + a{ul + ^ 2^0 — u^i — BY 

+ {1 — kl — akl)uQ + 2 {{ki — 1)A + ak 2 B — q) 

+ 2 mq {a{ui - U-iY + (1 - ki){u\ + u\) - ak 2 {u\ - vYY) . 

By requiring 

max {{ki — 1)^4 + ak 2 B) subjeet to 1 — kl — ak^ = 0, 

ki,k2 

we get 

ki = — , , /c2 = — . , 

YA^ + YA^ + aB^ 

and 

(fci — 1)^4 + ak2B = —A + Y A^ + aB'^. 

Now 21U beeomes 


21U = {ul + kiul + u\ — AY + Oi{ul + k2ul — u\ — BY 
+ 2 + YA^ + aB'^ — gj Mq 

+ 2mq (q!('Ui - U-iY + (1 - ki){ul + u\) - ak2{u\ - u\)) . (3.11) 
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3. For the last term, we elaim that 


a{ui — U-i) + (1 — ki){ui + u_i) — ak 2 {ui — u_i) > 0 for all ui, U-i G 
This quadratie form of mi, m_i ean be re-expressed as 

{a + 1 — ki — ak2)u\ + (a + 1 — fci + ak2)u\ — 2auiU-i. 


By using 0 < a < 1, we find that 


a + 1 — ki — ak2 = ( 1 — 


A 


+ a I 1 — 


B 


VA^ + aB^ J V 

a T 1 — ki T a/c 2 ^ 0, 

and the diseriminant of the quadratie form 


> 0 , 


(a + 1 — fci — ak 2 ){a -f 1 — fci -f a/c 2 ) — 

(1 + a)a^B^ 

(^2 + a52) [2^2 + a52 + 2A^/A^ + a#J 

This proves the elaim. 

4. For the term, we define the eonstant qi by 

gi := (-A + V/12 + aBA , 


> 0 


then the Mq term beeomes 


When q < qi, this term is positive. 


2 (gi - q) ul 


5. We have seen that for q < qi, all the terms of W are nonnegative and W = 0 only when 

“=(\/¥.o.\/¥)- 


6 . By plugging this eonstant state into the two eonstraints (1.2), (1.31, we get 


A = n, B = m. 


□ 
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3.3 Antiferromagnetic BEC {a > 0): qi < q < q 2 implies NS + 2C state 
Theorem 3.3. Suppose 0 < a < 1 and qi < q < q 2 . 

Ifm > 0, then the global minimizer of the variational problem (Pq) takes the form 


u = axf/ + bxQ\c/ 


(3.12) 


where 


SL = [\l ytind b = (0, + aB'^, 0 ), 

Tfl 

A = n + ir — l)q and B = —, 

r 

where r := |P|/|f2| satisfies 

2q^r^ + (2qn — q^) = 0 . 

Ifm = 0, then the global minimizer of the variational problem (Pq) has the form 

u = bxn where b=(0,A/n, 0). 


(3.13) 


(3.14) 


(3.15) 


Proof 1. Following the first two steps in Section [3^ we now choose A and B to cancel the 

e 

- A + VA^ + aB'^ = q. (3.16) 


Uq term in (3.111, i.e. we require 


Then 


2W{ui,Uo, M_i) = {ul + kiul + u\ — AY + a{ul + k2ul — u\ — Bf 


(3.17) 


+ 2ul [a{ui - u-ff + (1 - ki){u\ + u\) - ak 2 {u\ - u\)\ 


We have seen in Step 3 of Section 3.2 that the quadratic part in the last term is non¬ 


negative and it is zero only when ui = m_i = 0. Thus, from (3.17), hF > 0 and there are 
two roots for IF = 0: 


(a) Ko = 0 and from (3.171 «i, M-i satisfy the system 


u\ + u\ = A, 


ut — u 


1 


= B. 


(3.18) 


This is a 2C state. It takes the form 


u = 


-(y4 + P), 0, J-{A - B) 


= a. 
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(b) Ml = M_i = 0, but Mo 7^ 0. From (|3.17[), 


Mr 


= + a52. 


This is a NS state. It takes the form 


u = (0,^712+ a52,0) =b. 


In this ease, the minimizers of the variational problem 


inf / IF(mi, Mo, M_i) dx 

Jq 


take the form 


u — ^Xu + bxo\;7, 

where U is any measurable set in with relative size r := |[/|/|f2|. 

2. Our remaining task is to show the existenee of A, B and r for gi < g < q 2 from the eon- 
dition (3.16) and the two eonstraints ( |1.2| ) ( |1.3| ) and some natural inequality eonstraints. 
We list them below. 


—A + V A"^ + aB'^ = q 

(3.19) 

r A + {1 — r)\/ A^ + aB"^ = n 

(3.20) 

rB = m 

(3.21) 

0<r<l, 0<B<A, qi < q < q 2 - 

(3.22) 


The inequality 0 < i? < A is due to (eq:2C-AB). 

The first two equations give 

A = n + {r — l)q. 

The third equation leads to 

B = ^. 
r 

Substituting these two into the first equation, we obtain 

2q^r^ + (2qn — q^) — arri^ = 0 . 

We may rewrite it as the following dimensionless form 

+ (2q — — Tj = 0, (3.23) 

where q := q/n, rj = air? 0 < -g < 1. 
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Let us express the condition B < Am terms of r, q variables: 


1 ml 

7Tl ITl -L —- 

— = B<A = n + rq — q qr‘^ + (l — q)r ->0 q < h(r) := - 

r n 1 — r 


for 0 < r < 1. So our goal is to solve ( |3.23[ ) for x G [0,1] for given q G [qi, ^ 2 ] 
satisfying q < h{r). Here, q^ := qi/n, z = 1, 2. 


3. We rewrite (3.23) as a quadratic equation for q\ 


(2r^ — + 2r^g — rj = Q. 


(3.24) 


There are two branches of solutions for g: 

-r2 + ^r4 + z^(2r3 -r2) 


Qi{r) = 
Q2{r) = 


2r3 — 

—— \/r^ + ri{2r^ — r^) 
2r3 — 


(3.25) 


Since the turning point tq satisfying the equation + ri{2r — 1) = 0, we find tq = 
—Tj + \/rf + Tj. Furthermore, we have 0 < ro < \/2 — 1 < 1/2 because 0 < z] < 1. 

By direct calculation, we get that Qi{r) > 0 and decreasing on r > tq, and Q 2 {f') > 0 
and increasing on ro < r < 1/2. Furthermore, Q 2 {r) —?• cx) as r —)■ 1/2—. We plot the 
solution curve with n = 1, m = 0.2, a = 0.8 and the corresponding rj = 0.0032 in Figure 
[^. We also notice that h'{r) > 0 for 0 < r < 1 and h{l—) = 00 . Here, we have used 
m/n < 1. 


4. We claim that there is no solution for Q 2 {r) < h{r) with tq < r < 1/2. We calculate 

1 + \/l + (i - 


Q 2 {r) - h{r) = 


j, _ m 

n 


l-2r 




> 0 , 


(1 — 2r)(l — r) 

for ro < r < 1/2. Thus, there is no admissible solution on the branch q = Q 2 {r). 

5. Given q G [^ 1 ,^ 2 ], we look for r G [^o, 1] such that q = Qi{r) and q < h{r). We 
first notice that qi = Qi(l) from (3.25) and (3.8). On (ro, 1), the branch q = Qi{r) is 
strictly decreasing, and the function q = h{r) is strictly increasing. Thus, there exists a 
unique r 2 G (ro, 1) such that Qi{r 2 ) = h{r 2 ), because Qi{rQ) = Q 2 {ro) > h{rQ) and 
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Figure 1: n = 1, m = 0.2, a = 0.8, r] = 0.0032 


Qi(l)-Ml-) 
we get 


—oo. Indeed, at r 2 , we have A 


B. Then from (3.19), (3.20),(3.21), 


^ , -v/l + a — 1 / , /- , mx , . , 

Qi{r2) = - , (l + (a/1 + a - 1)—) = q2 = h{r2), 

V1 + a ^ ^ 

y/\ + d 

f2 = , - 

v^rT^-i + ^ 

For r G (r2, 1 ), we have Qi{r) < h{r) beeause Qi{-) is strietly deereasing and h{-) is 
strictly increasing. Now, we know Qi(l) = qi and Qi{r 2 ) = q 2 and Qi{-) is strictly 
decreasing in (r2), we get for every q G (gi, ^2), there exists a unique r G (r2, 1 ) such 
that q = Qi{r) and Qi{r) < h{r). 


6 . Lastly, we discuss the case of m = 0. We proceed the Step 1 and find that the global 


minimizer still has the form (3.12) and (3.13). We obtain x, A and B by solving (3.19), 


( |3.20| ) and ( |3.21| ) and get 

r = 0, A = n — q 


and 


B = 


(2n — q)q 


a 


This gives (|3.15[). 


□ 
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Remark 3.1. For antiferromagnetic EEC (a > 0), the 2C components Ui, m_i are suppressed 
by the value q. Consider the case that the total magnetization m > 0. In general, the larger 
value q is , the smaller values ui and m_i are. On the other hand, as q increases larger than qi 
the component uq begins to increase. When q is larger than q 2 , the component m_i vanishes and 
the ground state is NS + MS. 

3.4 Ferromagnetic BEC {a < 0): q < 0 implies MS + MS state 

Theorem 3.4. Suppose —1<q;< 0, m>0 and g < 0. Then the global minimizers of the 
constrained variational problem (Pq) tcike the form 

u = axt/ + b Xo\f/, 

where U is a measurable set of size 

\u\ = ^ (l + ™) IS^I 

and 

a=(yn, 0,0) and b = (0,0, ^/n)- 
Proof 1. First, we notice that 

‘IHj'p = + Uq + 'u_i) ~\~ Oi{u-^ — i) 3“ ‘^ouLqI^ui + M—i) + 2q(^Ui + u_-f) 

= (1 + a){ul + ul + - a (mq - 2|mi||m_i|)^ + 2q{ul + u^_i). 

Here we have used the following algebraic identity: 

2ul{\uf + \u.fY + K - - K + «o + «-i)' = -(«o - 2Kll«-i|)". 

2. Because of the constraint of total mass, we can convert our variational problem to an 

equivalent one by adding J^[—2{A{1 + a) + q){ul + Mq + + (1 + (S)A^] dx to the 

functional 2 J Htf dx. That is, the new energy density 

2W := 2Htf — (2^(1 + a) + 2q) {u\ + Mq + + (1 + 

= (1 + a){ul + Mq + u\ - AY - a (mq - 2|mi||m_i|)^ - 2qul (3.27) 

3. Since a < 0 and g < 0, every term of the function W is non-negative. The function 
equals to zero if and only if 

u = (v^, 0, 0) = a or u = (0, 0, v^) = b. 
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4. A measurable funetion u(a:) on f2 satisfies iy(u(a;)) dx = 0 if and only if there is a 
measurable setUcf2 sueh that 


u = ax{/ + b Xn\i7- 


5. By plugging sueh funetion u(a;) into the two eonstraints (1.2) and (1.31, we get 


A = n, \U\ = i (l + “) |n|. 


□ 


3.5 Ferromagnetic BEC (ct < 0): 0 < g implies 3C state 

Theorem 3.5. Suppose —l<a<0, m>0 and g > 0. Then the constrained variational 
problem (Pq) has a unique global minimizer 


U = 


where 



The value b is the unique root in (y/g^ + 2aqn, q) of the cubic equation 

b^ — (g^ + 2aqn)b + 2aq^m = 0. 


Proof. 1. As in the proof of the previous seetion, 

2Htf := {u\ + Mq + u\Y + a (2mq(|mi| + + {ul - + 2g(M^ + u\) 

= (1 + a){ul + ul + u^_iY - a [ul - 2|mi||m_i|)^ + 2g(M^ + 

Due to the two eonstraints, we ean add — 2a(M^ + Mq + m^i) and —2b[ul — u'^_f for 
some eonstants a, b to the above expression without ehanging the eonstrained variational 
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problem. We obtain 


2Htf — 2a{ul + Mq + u\) — 2b{ul — u\) 

= (1 + a){ul + Mq + M-i)^ - a (mq - 2|mi||m_i|)^ 

+ 2{q — a — b)ul — 2aul + 2(g — a + b)u‘ti 
= (1 + a){ul + Mq + - a (ul- 2|mi||m_i|)^ 

+ 2 [(g - a - bY^‘^\ui\ - {q - a + 

+ 4(g — a — by^^{q — a + — 2a|Mo|^ 


There will be two relations to determine a and b. First, we introduee the relation 

a = {q — a — — a + bY^^. 

This relation is equivalent to 


q b^ 


and leads to 


a = -, 

2 2g’ 


u ^ j, 

q — a — b = —- 0 = - > 0, 

^ 2 2q 2g - ’ 

, Q , (? + ^Y 

q — a + b= --l - h b = - > 0. 

^ 2 2q 2q - 


Sinee a should be non-negative from (3.28), we have 

q b'^ 

o = --->0 ^ | 6 |<,. 


(3.28) 


(3.29) 


(3.30) 


With this ehoiee of a, b satisfying the relations (3.29) and (3.30), we define 


2 iy = 2Htf — 2[a + (1 + a)n]{ul + Mq + u\) — 2b{u\ — vYY + (1 + a)n^ - 


= (1 + a){ul + Mq + u\ - riY - a (^ul - 2 |mi||m_i| + 


a 


+ - [iQ-b)\ui\ - {q + b)\u_i\Y 


(3.31) 


and the original eonstrained variational problem is equivalent to inf W (u(a;)) dx. 


2. For any given g > 0, fF of (3.31) has a unique minimizer (mi, mq, M-i). This leads to the 
following algebraic system for (mi, mq, M-i): 


2 1 2 1 2 
til + Wq + ^-1 ^ ^ 

(3.32) 

ul - 2|mi||m_i| = -- 
a 

(3.33) 

(g - b)\ui\ - (g + b)\u-i\ = 0 

(3.34) 
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For any fixed b, we solve this algebraic system for (mi, mq, m_i): 

\ui\ = 

\u-l\ = 


q + b 


1/2 

(3.35) 

[n + - 

) 

2q 

V 


q — b 

/ a ^ 

. 1/2 

(3.36) 

n + - 


2q 

V a/ 



l«o| = 


— b"^ + 6^ a] 

n + 


2g2 


2g2 a 


3. Our remaining task is find a relation to determine b. Since the constant state {ui, Mq, M-i) 
is the unique minimizer of W, we apply the constraint of total magnetization to this 
constant state and find 


uf — u\ = m. 


(3.37) 


From (3.37 ),(3.351 and (3.36), we obtain 


a 

n -\— 
a 


q + b 


q — b 


= m. 


Plugging p.29 1 into this equation, we obtain 

1 


or equivalently 


where 


b [n + ^^( 7 ^ ] = "^7) 


gib) = 0 , 


gib) := b^ — (g^ + 2aqn)b + 2aq^m = 0. 

This is the equation to determine b. In addition, there are other natural constraints that b 
should satisfy. In fact, subtracting p.33|) from p.32|), we obtain 


n-\ -= (|mi| + > 0 . 

a 


Combing this with (3.29) yields 


6^ > g^ + 2aqn. 


(3.38) 


(3.39) 


On the other hand, we have |6| < g from (3.30). Therefore, b must lie in the interval 

[^g2 + 2aqn, g]. 

4. We claim that g{b) = 0 has a unique root in [y/g^ + 4agn, g]. Because we have 

7 ( 3 / 7 ^ + 2aqn) = iq^ + 2agn)^^^ — (g^ + 2aqn)-s/q'^ + 2aqn + 2aq^m = 2aq^m < 0, 
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g{q) = — {q^ + 2aqn)q + 2aq^m = 2aq{m — n) > 0 


and 


g{b) = 36^ — + 2anq) = 2h^ + 6^ — (g^ + 2anq) > 0, 

The function g is strictly monotone on the interval [y/g^ + 2aqn, q]. Therefore, it has a 
unique root between on [y/g^ + 2aqn, g]. 

5. From the above discussion, we conclude there exists a unique 3C state when < 0 and 

g > 0. 

□ 


Remark 3.2. Suppose a < 0 and g = 0. In this case, following step 1 of Section |3.5| with 
a = 6 = 0, we set 


2W ■=2Htf — 2(1 + Q.)n(ul + Uq + mIi) + (1 + 

= (1 + a){ul + Mq + -nY - a (mq - 2|mi||m_i|)^ . 

The constrained variational problem (Pq) is equivalent to the variational problem 

inf / W{ui,Uo,u_i) dx. 

Jn 

Its constrained minimizer u = (mi, mq, m_i) satisfies 

u\{x) + Uq{x) + v?_i{x) = n 

- 2|Mi(a;)||M_i(x)| = 0 


for almost all X G kl. This gives |mi(x)| + |m_i(x)| = ^/n. Let us cdX\u\{x) — v?_i{x) = m(x). 
From the conservation of total magnetization, we should require / m(x) dx = M. We have 


Thus, we get 


(|mi(x)| - |m_i(x)|)(|mi(x) + |m-i(x)|) = m{x). 


n + mix) n — mix) 

“i(^) = —— ,u-i{x) = 


2\/n ’ 2\/n 

Since m{x) can be any arbitrary bounded measurable function with |m(x)| < n, there are 
infinite many Thomas-Fermi solutions in this case. 


3.6 BEC with {a = 0) 

Theorem 3.6. Suppose a = 0. Then the global minimizer of the constrained variational prob¬ 
lem (Pq) in a finite domain Vi is in either one of the following cases: 
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(i) Ifq = 0, then a minimizer takes the form u(a;) = (Mi(a;), Mo(a;), M_i(a;)) such that 

Ui(x) + Uq(x) + u^j^(x) = n for almost all x G kl with ul — u^_i dx = M. (3.40) 

Jq 

(ii) Ifq > 0, then a minimizer takes the form u(a;) = {ui{x), Uo{x), 0) such that 


ul{x) + ul{x) = n for almost all x G kl with / u\dx = M. (3.41) 

Jn 

(Hi) Ifq < 0, then a minimizer takes the form u(a:) = (ui(x), 0,M_i(a:)) such that 

ul{x) + vf_i{x) = n for almost all x G kl with / u\ — u^_^dx = M. (3.42) 

Jn 

Proof 1. When a = 0 and g = 0, we have Htf = + “o + We set 

W^(u) = - n{ul + 77 ^ + u\) + - = -(ul + ul + u\ - nf. 

Then the eonstrained variational problem (inf W(ui, uq, M-i) dx) is equivalent to the 

original one and its minimum is eharaeterized by ( |3.40| ). 

2. When a = 0 and q > 0, we follow the proof of Seetion [TTj In this ease, q 2 = 0 and W 
of ( |3.6| ) beeomes 

21V = (u^ + Mq + — Akf + Aquf-^. 

Beeause ly > 0 and the eonstrained minimization ean only oeeur when J W (u(x)) dx = 
0. This implies iy(u(a;)) = 0 for almost all x G kl. This forees U-i{x) = 0 and 
ul{x) + ul{x) + vj_i{x) = A. From the total mass eonstraint, we have to ehoose A = n. 

3. When a = 0 and g < 0, we follow the proof of Seetion [T4} In this ease, the normalized 
energy density W of (3.27) beeomes 

2W = {ul + Mq + u\ — AY — 2quQ. 

Similar to the argument of the previous step, we find uo{x) = 0 and A = n. 

□ 


4 P-convergence 

4.1 Interfacial and boundary energy functional 


The Thomas-Fermi solutions found in the last seetion are not unique in general. In faet, the 


pure states (Seetions 3.2, 3.5) are unique, while the mixed states (Seetions |3.1[ [X3 3.4) are 
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not unique. In the ease of mixed state, whieh has the form: u(a;) = axc/(a;) + bxn\;7(a:), 
only the ratio |?7|/|f2| is determined, but the measurable set U ean be arbitrary. It has been 
pointed out by Gurtin that interfaees are allowed to form without ehanging the bulk energy 
Htf(u(x)) dx [l26l[27]l . To seleet a physieal solution, we adopt the T-eonvergenee theory, 
whieh introduees an interfaeial energy funetional to penalize the formation of interfaees. This 
interfaeial energy funetional is the T-limit of the next-order expansion of the energy funetional 
£’e[u] as e —)■ 0. To be preeise, let us reeall that 

£’e[u] := / e^l Vup-fdx. 

Jn 

We write the domain of to be 

A:={ue n {L^{n,R+)f\Af[u] = N,M[u\ = M}. 


We expeet 


liminf £^e[ul = inf £^o[ul, 

e-S>0 


and thus look for the minimizers of So (i.e. the Thomas-Fermi solutions). Let us eall them 

Aq = arg min{£^o[u]| u G ]R+))^, A/'[u] = iV, A^[u] = M}, 

and the eorresponding minimal energy Eq. For mixed states, the set Aq, whieh is not a singleton, 
ean also be expressed as 


yio = {u = SiXu + bxo\c/| r = |t/|/|fl|, a, b are given in (|X3|), p.l4| ) or p.26D |. (4.1) 

We then define the next order energy funetional : dl —)■ M to be 

e.iu] = 

From the previous seetion, this funetional has the form 

Q^[u] := f e\Vu\‘^ +-W{u) dx, (4.2) 

Jn ^ 

where W is given in (3.4), (3.17), (3.27) whieh has the properties: iy(u) > 0 and lF(u) = 0 if 
and only if u = a or b. We expeet that 

lim inf [u] = inf Qq [u] . 

e—>'0 uG^l uG^Iq 

Here, the funetional Qq is so-ealled the F-limit of where the preeise definition will be given 
in Theorem |4.1 1 We will prove that : TIq —)■ is given by: 

2g{a, b) PerQ(u = a) for u = axu + ^Xn\u e A n {BV (12; M+))^ , 

+2g{0,a)'H‘^{{x G dO. : u(x) = a}) 

+2g{0,h)'H‘^{{x G dQ : u(x) = b}) 

-|-oo otherwise. 


^o[u] = < 
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where 


u) = inf 



: 7 : [0,1] —)■ Lipehitz eontinuous, 7 ( 0 ) 


V, 7(1) = u 


represents the minimal energy required to go from a eonstant state v to another eonstant state 
u. The notation Pern(u = a) is the perimeter of the set {u = a} in ff, and represents the 
two-dimensional Hausdorff measure. 

The intuition why Qq eontains an interfaeial energy ean be explained as the follows. It 
is expeeted that the minimizer of has a sharp transition from state a to state b aeross 
the interfaee dU fl fl, but have no variation up to order e along the tangential direction of the 
interface. The layer thickness should be of order e so that the kinetic energy e / | Vupand the 
bulk energy ^ f W(u) have the same order of magnitude. The minimal energy occurs only 
when these two energies are balanced, that is 



By the co-area formula, the energy contributed by the internal interface is roughly 2g(a, b)PerQ(u = 
a). This is the interfaeial energy. Similar argument can also explain the appearance of the 
boundary layer energy in Qq- 

Finally, the physical solution is selected by 


inf ^o[u]. 


ue.Ao 


This minimization problem is a geometric problem and can be solved by standard direct method 
in calculus of variations. 

4.2 Main Theorems 

We list our main theorems below. Although their proofs are mainly followed by the standard 
procedure of T-convergence arguments in [[55ll56ll . the quadratic constraints (i.e. A/^[u] = N, 
Ad[u] = M) in our present study require some modifications. We put these proofs in the next 
section for completeness. 

Theorem 4.1. The sequence {Q^} T -converges to Qo in -topology. This means the follows: 
1. (Lower semi-continuity) For any sequence {ue} C A converging to some uq G Aq in 


we have 


^o[uo] < liminf^Ju,]. 
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2. (Recovery sequence) For any vq G Aq, there exists a sequence {v^} C A converging to 
vq in such that 

^o[vo] = lim^Jv,]. 

e-s>0 

Such Qo is called the T -limit of 

Theorem 4.2. Suppose that {u^} is a family in A with an uniformly bounded energy, that is 

Ge[n,] < Co (4.3) 

for some positive constant Co- Then there exists a subsequence {u^^} converges to some uq G 
Aq in as j —)■ oo. 

Theorem 4.3. Suppose are minimizers of the variational problem 


inf Gfu] 


and {u^j } converges uq G .Aq in (Lf(FL))^ for some subsequence Cj —)■ 0. Then uq solves the 
variational problem 

inf ^o[u]. 


Proof Let vq G Aq. There exists a sequence G A, such that ^o[vo] 
have 

^o[vo] = lim^Jv,] > liminf^Ju,)] > ^o[uo]. 

e—^-O e—^0 

This shows that uq minimizes the functional Go- 


lime^o^e[ve]. We 


□ 


According to Theorem 4.3 and the discussion in sections |3.1[ |3.2[ |3.3[ 3.4 and 3.5 


we 


conclude that the asymptotic behaviors of the ground states of the Spin-1 BEC systems are 
characterized by their corresponding Thomas-Fermi solutions satisfying the minimal interface 
criterion. 

By the work of ll3^ . it is also possible to construct local minimizers of the perturbed varia¬ 
tional problem from isolated local minimizers of the limiting one. 


Theorem 4.4. Suppose Uq G Ao is an isolated -local minimizer of the functional Go, i-e. 
there exists e > 0 such that 


^o(v) > ^o(uo) for all v G A,o and 0 < ||v — uo||l 2 < e. 

Then there exists a sequence {ue} C A such that each is a local minimizer of the functional 
Ge and Ug —>■ Uq in ai e —)■ 0. 
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d{du nn) = d{du n dn) 


Figure 2: The geometry of the domain and the interfaee dU fl 


The existenee of a minimizer for the limiting problem 

inf ^o[u] 

ue^lo 

is obtained through the standard direet method in the ealeulus of variations. By the straight 
forward ealeulation of the first variation, we get the following neeessary eondition for the inter¬ 
faee. 

Theorem 4.5. Let flbe a bounded domain with a -boundary dVL and uq be a critical point of 
Qq such that dU fl fl is of class with mean curvature H : dU fl fl —)■ M. The corresponding 
Euler-Lagrange equation o/min^o[u] is 

J H{x) = Const. for X G dU fl fl, 

I g{sL, ■ t) + 5f(0, a) — g{0, b) = 0 for x G d{dU fl fl) 

where v : d{dU fl fl) —)■ is an outward unit tangential vector to the interface dU fl fl and 

normal to d{dU fl fl); t : d{dU fl fl) —)■ is the outward unit tangential vector to dU fl dfl 

and normal to d{dU fl dfl). 
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Proof. The calculation the first variation of Qq under the fixed volume constraint is similar to 
mm. We omit it here. □ 

Remark 4.1. The second equation of the Euler-Lagrange equation is called the Young’s relation, 
which appears in natural process of wetting: 

g{aL, b) cos 6* + 5 f( 0 , a) — 5 f( 0 , b) = 0. 

Here, the contact angle between the boundary and the interface is denoted by 9, see Figure 
Indeed, it is a balance law between surface tensions of three different interfaces on the boundary. 
Thus, the Euler-Lagrange equation mentioned above could also be considered as a equation for 
a quantum-like wetting process. 

4.3 Preliminary lemmas 

The function IE : —)■ M constructed in those sections has the following properties: 

1. IE is a C^-nonnegative function with the following symmetry property: 

IE (mi, Mo, m_i) = IE (-Ml, Mo, m_i) = IE(mi, -mo, m_i) = IE(mi, mo, -m_i) (4.4) 

and 

IE(u) = 0 in if and only if u = {a, b}. 

2. There exist 5 > 0 and C > 0 such that 

C|u — a|^ < IE(u) < —|u — a|^ for |u — a| <6 

o 

and 

Clu-bp < IE(u) < ^|u-b|2 for|u-b|<(5 

(_y 

3. There exist two positive values C and R such that 

C|u|2 < IE(u) for |u| > R. (4.7) 

We shall assume these properties of IE in the discussion below. 

We quote several useful lemmas from |j3l |43l |551 HH which will be used in the proof of our 
r-convergence result. 


(4.5) 

(4.6) 
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Lemma 4.6 (See B3ll55]| ). Let f2 be an open bounded subset in M” with Lipschitz-continuous 
boundary. Let A be an open subset in M” with compact -boundary dA such that H‘^~^{dA fl 
dfl) = 0. Define the signed distance function to dA, d \VL by 


dA{x) 


dist{x, dA) X G fi\ 74 , 
—dist{x, dA) X G A n 


Then for some s > d is function in {|(i(x)| < s} with 


(4.8) 


IV^aI = 1. 


Furthermore, 

\imn^-\{x eVt- dAix) = s}) =H'^-\dAnVt). (4.9) 

s^O 

Lemma 4.7 (See (Si ESI EH). For any G M+, there exists a curve : [0,1] —)• 

with 7 ( 0 ) = and 7 ( 1 ) = which minimizes the degenerate geodesic problem 

■= inf y/Wi-tir)) | 7 '(r)| dr : j e Lip([ 0 , 1 ]; R^), 7 ( 0 ) = 7 ( 1 ) = . 


Define (^^(,(^ 1 ) = Then the function —)■ M+ is Lipchitz-continuous with the 


property 

^ for almost all $,1 G M+. 

Furthermore, ifu G H^{il;M.^) fl L°°(f2; Mij_), then o u G iy^’^(f2; Mi^) and 



o u) I (ix < 



\/iy(u(x))|Vu(x)| dx, 


where | Vu| is the 2-norm o/Vu, that is 


Vu 



(4.10) 


Proof 1. This geodesie problem eorresponds to a Riemannian metrie ds^ = Wds\^^ on 
where dsEuc is the Euelidean matrie. However, there are two diffieulties here. The 
first one is that this is a eonstrained variational problem, namely all eomponents of the 
geodesie eurve 7 (-) should be non-negative. To resolve this diffieulty, we extend this 
variational problem to the entire by taking the advantage of the symmetry property of 
the funetion W. Thus we eonsider an equivalent degenerate geodesie problem on 

inf \/W{'p{t)) \i{T)\dT : 7 G Lip{[0,l];R^), 7 ( 0 ) = ^ 0 , 7(1) = ^i| , 

(4.11) 
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For any ^q, e if we have found a Lipschitz geodesic 7 (-) connecting to in 
then using the symmetry property of W and reflection, we can always find a represen¬ 
tative Lipschitz-continuous curve 7 : [0,1] —)■ which solves the constrained geodesic 

problem. 


2. The second difflculy is the degeneracy of W{u) at a and b. Thus, the direct method in 
the calculus of variations cannot be applied straightforwardly. Therefore, we consider the 
regularized problem: 


inf|y : 7 e Lip([ 0 ,l];M^), 7 ( 0 ) =^o, 7 ( 1 ) 

(4.12) 

where 5 > 0. The regularized Riemannian metric dsj = {W + is conformal to 

the Euclidean metric on the plane. Its minimizer 7 ^ exists uniquely by the direct method 
in the calculus of variations. A uniform bound on the Euclidean arclength of 7 ^ can be 
derived by carefully analyzing the curve 75 in two different regions: one region is away 
from the zeros of W and the other is the region near the zeros of W. Because the value 
Is is invariant under re-parametrization of the curve 7 ^, this allows us to choose a new 
parametrization 7 ^ : [0,1] —)■ with a constant speed, that is I 75 I = ^s. By the Arzela- 
Ascoli compactness theorem, there exists a subsequence { 7 ^^} and a curve 7 such that 
7 r converges to 7 in C([0,1], M^). It is easy to see that 7 solves the problem (|4.11[). 


3. It is observed that g : x —)■ M is a metric in and satisfies the triangle inequality: 




That is, 

V?4o(^i) ^ ^ ( 74 ;, 4 i(^))1 7474 ,(^)(4.13) 

Erom this, we get that is locally Eipchitz continuous and 

|Vv>£,(«i)l < (4.14) 

As we choose moving along the geodesic 74 ^, 4 ^ to approach then the inequality in 
( |4.13| ) becomes equality, and we get 

(4.15) 

Notice that the above arguments hold for all in the interior of But the inequality 
(4.14) and equality (4.15) can be extended to the boundary of 
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4. Applying the Cauchy inequality and from (|4.14[), we get 


/ ou)|cix = / |(V(^^g(u)) V^u| dx 

Jq Jq 

< / I V(^ 4 (j(u)| iV^ul dx < / a/1C(u(x))| Vu(x)| dx 

Ju Ju 

foru G By the density theorem, this inequality also holds for u G 

□ 

The following lemma is used to construct the one-dimensional profile of the internal layer. 

Lemma 4.8 (See [l55ll56l l. Given a, b G Then there exists a Lipchitz-continuous function 
T] : (— 00 , 00 ) —)• whose trajectory is the geodesic with the metric ds^ = Wds%^^ connecting 
a to b. The function r]{-) solves the variational problem 


inf / 

■n{-oo)=a J 
T7(+oo)=b 


r/'f + W{r]) ds 


with minimal value 2g{sL, b) and 

\r]it) — a| < as t ^ — 00 , 

\r]{t) — b| < as t ^ +cx). 


(4.16) 

(4.17) 


Proof Let 7 (-) be the geodesic curve in Lemma [4?^ connecting a to b. Let us parametrize it 
by G [0,1] such that 7 ( 0 ) = a, 7 ( 1 ) = b, 7 (-) is Lipschitz continuous and | 7 '(/ 9 )| > c > 0 
for some constant c and for all [3 G [0,1]. Now, let us consider the ODE: 

dp _ \/W(7(/3)) 

dtlO - |Y(/3)| , 

/3(0) = i. 


The right-hand side of this ODE, /(/3) := /\jf(3)\, is Eipchiptz continuous on 

[0,1], thus we have local existence and uniqueness of the solution. Note that /3 = 0,1 are 
the only two zeros of /, thus, from uniqueness, the solution of this ODE stays between 0 and 
1, as long as it exists. Thereby it exists globally. Eurthermore, j3{t) —)■ 0 (resp. j3{t) —)■ 1) 
exponentially fast as f —)■ —00 (resp. t —>■ 00 ) because /'(O) (resp. /'(I) 7 ^ 0), which is, in turn, 
due to ( |4.5[ ) (resp. ( |4.6[ )). 

Now let us define r]{t) := 'y{/3{t)). We have 


1 ^ 


'(i)i = l7'(^^(^))^^'(^)l = h'wmptt) = = s/wMf). 
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Let ff G Lip{{—oo, oo), be any eurve eonneeting a to b. Using the Cauehy inequality, the 
faet that 7 is geodesie and \r]'\ = W{r]),w& get 



+ W{fi) dt 


> 2 



\/W{fi)\f}'\ dt >2 



\/W(^\r}'\ dt 



+ W{r}) dt 


Thus, 77 solves the variational problem. 


□ 


Similarly, we also have the lemma for the eonstruction of the one-dimensional profile of the 
boundary layer due to the differenee between the Diriehlet boundary condition u = 0 on dU 
and u = a, b in U. 


Lemma 4.9. There exist two Lipchitz-continuous functions : (0, 00 ) —)■ and : 

(0, 00 ) —>■ whose trajectories are the geodesics with the metric ds^ = Wds\^^ connecting 

0 to a and to b, respectively. They solve the variational problems 


inf / + W{r})dt 

»?(o)=o Jo 

T7(+oo)=a 

and 

poo 

inf / Ir/'P + lUfr/) dt, 

^(o)=o Jo 

T 7 {+oo)=b 

with minimal values 2g{Q, a) and 2g{Q, b), respectively. Furthermore, limi_)._|_oo 'naif) = ^ 
linit_^+oo T/b (f) = b are being attained at exponential rates. 


Proof. The proof is similar to the proof of Lemma [4^ □ 


4.4 Proofs of the Main Theorems 


Proof of Theorem 14. II 


Proof of Lower semi-continuity. 


1. Let us extend and uo trivially to a larger bounded smooth domain U' such that U CC U'. 
That is. 


u,(a:) 


Ug(a;) forxeU, 

0 for a: G 


and Uo(x) 


uo(a;) forxGU, 

0 for X G 


We have —)■ uq in This together with the fact that the two functions and po are 

Lipschitz continuous ( Lemma [477]) lead to 


o w)(a:) iPa o Uo)(a:) 


0 X eU 

5 f(a, b) X E fl\U inL^(U') 

gisL, 0) X E Q'\fl 
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and 


g{0, a) X ^ U 

((^0 o u,)(a;)-)■ ((^0 o uo)(x) = <( 5f(0,b) x E ^l\U mL'^{W). 

0 X E f2'\f2 

2. Let and be two bounded sets defined by 

A^ = {x e : dist(a;, dO.) >5} and A2 = {x G : dist(a;, < 5} U (f2'\f2). 


By using the inequality of arithmetic and geometric means and the lower semicontinuity of the 
BV-norm under L^-convergence, we have 

liminf > 2 liminf / \/W{\i^{x)) |Vu e(a;)| dx 


e^O 


e-S>0 


= 2 lim inf 

e-5-O 




v/lL(fi,(a;))|Vri,(a;)|da;+ / ^/W{n,{x)) |Va e(a;)| dx 


lAl 


>2 lim inf j |Vue(a;)| dx + 2 liminf j ^lL(ue(x)) |VUe(x)| dx 


lAi 


JaI 

After applying the inequality ( |4.10| ) to the function (p^ o on the domain Af and the function 
ifo ° Ue on the domain A 2 , we find that 


lim inf ^JuJ > 2 lim inf / iVf^^a o Ue)| dx + 2 liminf / iVfv^oouAldx 

e-5-O e-!>0 J^s £-5-0 J^s 


> 2 


lAf 


|V(99aOUo)| + / |V((^oOUo)| . 


lAi 


Taking d —)■ 0, we get 


lim inf > 2g{aL, b) Pern(u = a) + 2g{0, a) 'H^({x E dQ : u = a}) 

e—^0 

+ 2g{0, b) H'^{{x E d^l : u = b}). 


□ 


Proof of Recovery sequence. 

1. We will construct recovery sequence {vg} C A for vq G Aq (4.1) of the form 


Vo = 


if X G V, 

ifx G n\v, 


(4.18) 


and satisfies A/'[vo] = A^, Ad[vo] = M. We discuss the case when vq G Aq H 

first. The case when vq G Aq \ {BV will be discussed in item 6 below. In the first case, 

since vq G {BV{il))^, the set V has finite perimeter in 12. We may assume that the interface 
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r = dV n f2 is smooth because a set of finite perimeter can be approximated by a sequence of 
sets with smooth boundary, see ll55ll . 

2. The recovery sequence to be constructed will have the form 

v,(a:) = w,(a:) + ai(e) ip{x) + a 2 {e) 


Here, is a layer solution; ^p{x), '0(x) are smooth functions supported on V and il\V, 
respectively; the coefficients ai(e), 0 : 2 ( 6 ) are 0{e). The terms oi(e) (p{x) + 02 (e) jJj{x) are 
designed so that the conservation constraints = N and Ad[ve] = M are satisfied. Detail 

conditions for (fi{x) and ipix) will be given later. The layer solution involves an internal 
layer function rj and an auxiliary boundary layer function 77 . The internal layer function t) 
is defined in Lemma 4.8[ which connects a to b. The auxiliary boundary layer function fj : 
{— 00 , 00 ) —)■ is defined by 


fi{t) = 



if f < 0 , 
if f > 0 , 


where 77 ^ and are defined in Lemma 4.9 They are the boundary layer functions connecting 
0 to a and to b, respectively. Note that the function 77 (f) on the negative (resp. positive) half 
line {t < 0 } (resp. {t > 0 }) is the profile of the boundary layer approaching to the value a 


(resp. b). Let us also define the signed distance function d = dy associated with dV by (4.8) 
and an auxiliary distance function associated with dVt by 


db{x) = 


-dist(a:, dVt) if x G V, 


dist(x, dVt) if X e D \ L. 

We choose a cut-off function C,, which is a smooth function such that 0 <c< 1 and 


m = 


0 


if |f| < 1 , 
if |f| > 2 . 


Let us set 7 = |. Finally, we define the layer solution by 


wJx = 


1-C 

+ c 


dh{x) 

e')' 

db{x) 


C 

i-C 


d{x] 

d{x 

e'y 






db{x) 


1-C 


d{x) 


Vo (a:) 


3. We claim that the sequence {w^} converges to vq in (L^(D))^ and 

A/’[wg] = iV + rig, Ad[we] = M -\-m^ 


(4.19) 


33 





















Figure 3: The domain represents the domain ad the region on the left side of the solid line 
dV n represents the domain V. The distance between each type of line is e"'. The fc-region 
represents the sub-domain flfc. 


with = 0 (e) and = 0 (e). 

We partition the domain Q into sub-domains Qk (See Figure [^: 

= {a: G : \d(x)\ < e^, |db(a:)| > 2e'’'}, 

= { 2 : £ (2 : < \d(x)\ < 2e'^, \db(x)\ > 2e^}, 

0,3 = {x e Q : \d(x)\ > 2e'^, \db(x)\ < e'’'}, 

= {x G fl : |d(x)| > 2 e"',e'^ < \db(x)\ < 2e'^}, 

fls = {x G fl : |d(x)| > 2e'^, \db(x)\ > 2e'^}, 

fie = ^ : \d(x)\ < 2e"', \db(x)\ < 2e'^} 

and estimate the following integral term by term on each subdomain (1^: 

/ |w, -vo|dx = ^/ |w, -vo|dx. 

J Q . fc=l ^ 


34 






We calculate 


/ I We — \o\dx < / 

' {xeniUn 2 : 0<(i(a:)<2eT'} i{0<d(x)<2eT'} 




< 


' { 0 <d{x)< 2 e'y} 
r2e'< 


e-y 




e 

dx 


dx 


r]{-)—h {d{x) = s} ds 


< e ( max T-L^ldix) = s} 

0<s<2eT' 


r- 2 eT'- 


|r/(t) -b| 


because the exponential decay estimates (4.16), (4.17) and (4.9). Similarly, we have 


'{3;eOiUr22: - 2 e^<d(x)< 0 } 


We—vo| dx < e i max 'H‘^{d{x) = s} 


-2eT'<s<0 / _/_2e7-l 


|r/(t) -a| c 


'{3;g03U04: — 2 e'> <di,{x)< 0 } 


|we—vo| dx < e {dh{x) = s}^ J l^a(^) ~ ^ 


-2eT'- 


7{3;eO3Ur24:0<d(a:)<2e')'} 

We calculate 


and 


Finally, we calculate 


|we—vo| dx < e ^^max 'H^{d{x) = s}^ J |^7b(^) ~ b| dt 


I We — Vo I dx = 0 


'Os 


/ |we — vo| dx < max(|we| + IvoDl^lel = 0{e^'^). 
iQe 


|WeP — IvqP dx 


(|We| + |vo|)(|We| - |wo|) dx 


In 


Similarly, we have 


< max(||we|| + ||vo||) / |we - vq] dx = 0 (e). 

® Jn 


(ime,ir - k.,-ir) - (ko,ip - ko,-ir) dx = o(e). 


4. We rewrite ^e[we] and estimate them term by term: 


^JWel = 


f e|VWeP + -hF(We) dx = ^ [ e|VWeP + -PF(We) dx. 

Jn ^ k=i J^^ ^ 


dt = 0(e) 

It = 0(e), 

-I dt = 0(e), 
= 0(e). 
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Using the coarea formula, ( |4.9| ) and Lemma [4T8| we calculate the first term 

1 , 


lim / e|Vwe| +-lU(we) dx 


e-S-O 




= lim - 
£—>■0 e 


= lim - 
e-!>0 e 


1 




, f d{x) 


+ W(r] 


d{x) 


dx 


= lim 

e-S>0 


r 

'-£7-1 


/ t S\ 


^ (“) 

+ 1U (^ (-)) 

Ve/ 

V Ve/ / 


G Ui : d{x) = s} ds 


(|77'(f)p -f W{ri{t))) 'H‘^{x e Ui : d{x) = et} dt 


= 2g(a, b) lim max 'M?{x G U : d(x) = s} 

c-S-O |s|<eT' 

< 25((a,b)PerQ(w = a). 


Similarly, for boundary layer, we also have 

1 , 


lim / elVwel -f-lU(w£) dx 


£^0 


'Os 


= lim - 
£—>■0 e 


O 3 




, f db{x) 


+ W\f, 


db{x) 


dx 


= lim- f 

e-!>0 e 


-eT 

7-1 


/s\ 

2 ,,, /s\\ 

^ (“) 

+ (^(-)) 

Ve/ 

V Ve// 


G U 3 : db{x) = s} ds 


= lim 

e-S>0 


= lim 

e-S>0 


-e7-l 


(|r/'(f)p -f W{ri{t))) "H^jx G U3 : db{x) = et} dt 


e U3 : db{x) = -et}dt 

(iVbWP + ^(^b(^))) e U3 : db{x) = etjdt 


. 7-1 


< 2gf(0,a) lim max T-L^tx G U : db(x) = —s| 

“ ^ ^ £^0 0<s<£'r 0 \ J ) 

+2g(0, b) lim max H^jx G U : db(x) = s} 

= 2g{Q, a) 'H‘^{{x G dVL : Vo(x) = a}) -f 2g{Q, b) "H^dx G dVt : Vo(x) = b}). 

Applying the Taylor expansion of the function W around a and b and using the exponential 


decay estimates ( |4.16[ ) and ( |4.17[ ) 

lim / elVwgp-f-fUfwe) dx = 0 , 

do, ^ 


lim / elVwgp-f-fUfwe) dx = 0 . 

d04 ^ 


36 























Since equals to a or b on fls, we have 


e|Vwe|^ + -iy(we) dx = 0 (e). 


'05 


Thus, we obtain 


e| Vw ,|2 + -iy(w,) dx = 0(e) + 0 (€^^-^) 


'fig 


limsup^Jw,] < ^o[vo]- 

e-S>0 

Combining the result of the lower semi-eontinutiy, we have 


lim^Jw,] = ^o[vo]- 

e^O 

5. Finally, we modify the layer funetion by adding some smooth function with compact 
support in order to fit the conservation eonstraints. We choose two smooth functions : fl —)■ 
and t/? : fl —)■ such that t/? satisfy the following eonditions: 

(i) The function ip has compact support in V and the function 1 /? has compact support in 


(ii) There exists 5 > 0 sueh that all the eomponents of the funetion (a + 'y(p(x)) and the 
funetion (b + 'j'ijj(x)) are all nonnegative on fl whenever I 7 I < 6 . 


(iii) The matrix 

Jyiaipi + aoipo + a-iip-i) dx + &o' 0 o + dx 

Jyiaiifi - a-iip-i) dx dx 

is invertible where a = (oi, Oq, a_i) and b = (foi, feoi ^-i)- 

For eaeh e small enough, we would like to find ai, 02 sueh that the funetion (w^ + aicp + 
Q!2'0) satisfies the eonstraints of total mass and total magnetization. That is 

\we,i + aipi + a2'ipi\^ + \wefl + aipo + a2'?/’oP + \we-i + aip-i + a2'rp-i\^ dx = 

Jq \we,i + aipi + a2'?/’ip - Itr'e-1 + + a2'?/’-iP dx = 


Because of (4.19), we obtain the system 


/i(q!1,Q!2) = «! / ((Pi+ (Po + (p_i)dx + 2 ai / (We,l(pi + W^flPo + We,l(p-l) dx 
Jn Jn 

+«2 / (V’l + V'o + V’-i) dx + 2a2 / (w^^i'ipi + w^^oi’o + dx = 

Jn Jn 

12(^1,02) =al / ((pi - dx + 2 ai / {we,ipi - We-iP-i) dx 
Jn Jn 

+Q!2 / (V’l ~ V'-i) dx + 2a2 / (we^iipi — We-i'ip-i) dx = —m^. 


-n. 
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It is observed that (ai, 02 ) = (0, 0) is a solution of the system 


f /i(ai,a2) = 0, 

\ /2(ttl,tt2) = 0, 

and the Jaeobian matrix 

f Ift(O’O) 2^ + + + bo^o + dx 

V IS(O’O) / \ fy(am - dx f^^y(bi^i - b-i^-i) dx 

is invertible. Aeeording to the Inverse function theorem, for e small enough the system 

f fi(o:i,a2) = ~ne, 

\ f2{ai,a2) = -me, 

is solvable. Furthermore, for each e small, the corresponding Q;i(e), 0 : 2 ( 6 ) are of order 0(e^'^). 
Define the recovering sequence {v^} by 

Veix) = We{x) + ai(e) cp{x) + 02(6) 'ljj{x). 

By the choice of the function v^, it will satisfy the constraints of total mass and total magneti¬ 
zation. Finally, we calculate 


e |V(we -f oi(e) (f + 02 (e) i/?)| + -W{we + oi(e) ip + 02 (e) 1 /^) dx 


e^O 


e^O 


/ 2 

e |V(we -1- oi(e) ip + 02 (e) i/’)! -f -W(we + oi(e) p + 02 (e) i/>) dx 




'05 


lim 


e |V(we -f oi(e) p + 02 (e) + -W{we + oi(e) p + 02 (e) 1 /^) dx 


Jn\n 5 e 

When e is small enough, the support of the function p and 1 /^ is contained in fls. We find that 


/ 2 

e |V(we + oi(e) p + 02 (e) i/’)! + -W(we -f oi(e) p + 02 (e) 1 /^) dx 




'05 


= lim / -ll^(we + oi(e) cp + 02 (e) i/^) dx 
< Clim [ -1 oi (e) + 02 (e) 1 /^ 1 ^ dx 


/05 


= 0 (e). 

Since Oi(e) + 02(e) 1 /^ = 0 on f 2 \f 25 for e small enough, we obtain 


lim^Jv,] = ^o[vo]- 
€^0 
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Thus, {vg} is a recovery sequence for the functional Qo about vq G D BV(Q). 

6 . Finally, we need to construct a recovery sequence {vj C A for the case when vq G 
Ao\BV(Q). This means that the set V := {vq = a} has an infinite perimeter. Our approach is 
to construct a sequence {v,;} of the form 

Ve(a;) = We(a;) + Q;i(e) cp(x) + 0:2(6) 't/y(x), 

then to choose IDqI; Here, the sequence {w^} C is taken to be 

We := Vq * 0 e, whcrc 0e IS a standard modifier; and cp : Q —)■ and ; O —)■ are two 
smooth functions chosen in a way that the matrix 

/ /^(v^Pi + vIIpo + dx + Vo'ipQ + dx \ 

V jM^^-v\p-i)dx j^{v\i)i-v^_^i)-i)dx ] 

is invertible. The coefficients Oj(e) are chosen so that Ve(a;) satisfies the conservation con¬ 
straints. It is clear that —> vq in (L^(f2))^. We thus get 

A/”[We] —)■ A/”[vo] = iV, A^[We] —)■ [vq] = M. 

By the implicit function theorem, we obtain that there exists some eo > 0 and two corresponding 
functions ai(e), 02 (e) for 0 < e < eo with Oj(e) —)■ 0 as e —)■ 0 such that the corresponding 
function satisfy 

Ar[v,] = iV, Af[v,] = M. 

From Oj(e) —)■ 0 as e —)■ 0 and w,. —)■ vq, we get —)■ vq in (L^(f2))^. Now set = 
(Ifi^l, iFoU^iil)- It is obvious that the sequence {v^} also converges to vq in (L^(f2))^ and 
satisfies the conservation constraints 


A/'[v,] = iV, Af[v,] = M. 


By the lower-semicontinuity of Qe that has been proved in the first part of this theorem, we have 


-fcx) = 2 g{a, b)Per(vo = a) < ^o(vo) < liminf ^,(v,). 


This implies that 


lim^£(ve) = -fcx) = ^o(vo)- 

This completes the proof of the F-convergence of the sequence {G^} to Go- 


□ 
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Proof of Theorem |4.2| Set Ri := max(i?, |a| +1, |b| +1) where R is given in ( |4.7[ ) and define 
a truneating funetion by 


By using (4.7) and the assumption (4.31, we find 



lUfT dx 




1 

< — 
- C 


{|Ue|>r?l} 


dx 




(4.20) 


as e —)■ 0. Therefore, we only need to show the precompaetness of the L°°-sequenee {ve} in 
(L^(f2))^. In order to aehieve this, we will apply the eompaetness result for Young measures: 
the L°°-boundedness of the sequenee {v^} implies there exist a subsequenee and a Young 
measure /i sueh that 


if / is a eontinuous funetion on then /(vg) ^x ^ J f{y) djjixiy)^ in weakly*. 

(4.21) 

In the light of ( |4.20[ ), we eonelude that the measure of the set {x G fl : |u£(x) | > i?i} tends to 
zero. We find 


lim / W (v,;) dx = lim 


e-5-O 


In 


e-5-0 


'{|u,|<r?i} 


IT(uJ dx + 


'{|Uf|>Kl} 


iy(0) dx 


< lim / W{u^)dx 


e-5-O 


< e ^e[Ue] —>■ 0. 


Because of (4.4), we know that 


fix = d{x) 6y=a + (1 — d{x)) 6y=b a.e. x E fl 


(4.22) 


where 0 < 9{x) < 1. 

On the other hand, from Lemma we know that is Lipchitz continuous and leads to 
the -boundedness of the sequence o v^}. Thus, there exist a subsequence o } and 

a Young measure Ux such that the subsequence {(^a o converges to a Young measure Ux. 


Because of (4.22), we can express the Young measure Ux as 


Ux = 0{x) Sy=^^fa) + (1 - d{x)) a.e. x eQ. 
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Now, we are going to show the sequence { 99 a ° is bounded in BV (f2). First, we estimate 

/ |(^a(u,(a;))| dx < / |(^a(v,(x))| dx + / |(^a(Us(a;)) - (^a(v,(x))| dx 

Ju Jq Ju 

< |fi|||</9a O vJloo + / \(p'^\\u^ dx 

Jq 

and find (fa o is bounded in since the function ipa{x) is Lipchitz continuous on and 
||ue — Ve ||£,2 —)■ 0 as e —)■ 0. Next, we estimate 

I |V((^a o uj(x))| dx < I ^iy(u,(x)) IVu^l dx 
Jn Jn 

< - f elVUfp + -hF(uJ dx 

2 Jo e 

< < Y < +00- 

Therefore, the sequence {(^aOUg} is bounded in BV-norm. There exists a subsequence (without 
abusing the notation, we still denote it by the same sequence) and a function h G such 

that 

(^a ° Ug ^ inL^(f2). 

Since the function (^a is Lipchitz continuous on and ||ue — V£||j ;^2 — 0 as e —)• 0, we have 

</9a o Ve —)■ h inL^(f2). 

In consequence, the Young measure h>x associated with (^a o is just a point mass a.e. x G fl, 
that is 

Vx ^h{x) a.e. X G fl 

and the function 6 {x) = xu for some measurable set (/ C fl. Thus, the corresponding Young 
measure could be represented by the function 

Uo := axf/ + bxo\;7- 

The definition of the convergence in Young measures to a function uq give us that 

Ve —)■ Uq weakly in U’iyt) for 1 < p < cx) 
and 

llVelliP —)■ ||uo||lp for 1 < P < CX). 

It follows that V,; converges to uq strongly in (L^(f2))^ for 1 < p < 00 . From ( |4.20[ ), we have u,. 
converges to uq in (L^(f2))^. Because Uq also satisfies the conservation constraints A/'[uo] = N 
and A^[uo] = M., we can conclude uq G Jlo- This completes the proof. 

□ 
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